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•  An  electrochemistry-based  impedance  model  for  lithium-ion  batteries  is  proposed. 

•  An  optimization  method  is  used  to  identify  the  model  parameters  from  EIS  data. 

•  The  accuracy  of  the  proposed  model  is  compared  with  that  of  the  Randles  model. 

•  Some  parameters  of  the  proposed  model  show  clear  trend  with  aging  cycles. 
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Accurate  models  of  lithium-ion  batteries  are  important  for  analyzing  and  predicting  battery  dynamics 
and  aging.  This  paper  presents  an  electrochemistry-based  impedance  model  for  lithium-ion  batteries  to 
better  understand  the  relationship  between  battery  internal  dynamics  and  external  measurement.  The 
proposed  impedance  model  is  a  modified  single  particle  model  which  balances  between  simplicity  and 
accuracy.  The  model  includes  electrochemical  impedance  due  to  charge-transfer  reaction,  diffusion  dy¬ 
namics  in  the  electrodes,  effects  of  ion  concentration,  capacitance  dispersion  in  the  double  layer,  and 
anode  insulating  film  growth,  etc.  The  impedance  tests  for  model  validation  were  performed  on  two 
lithium-ion  cells  at  ambient  temperature  and  at  different  SOC  levels.  A  particle  swarm  optimization 
method  is  employed  to  identify  model  parameters.  The  model  accuracy  under  different  conditions  is 
compared  with  that  of  conventional  Randles  model  and  the  parameter  variations  at  different  stage  of  the 
aging  process  are  studied. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  are  promising  energy  storage  devices  for 
portable  electronics,  electrified  vehicles  and  renewable  energy 
storage  [12].  The  prediction  of  battery  performances  is  important, 
especially  for  highly  dynamic  applications,  e.g.  electrified  vehicles 
and  grid  storage  systems  [3,4].  In  the  past  decades  considerable 
efforts  have  been  made  to  analyze,  model  and  control  battery  be¬ 
haviors,  including  open-circuit  potential,  charging/discharging 
dynamics,  capacity  loss,  etc.  [5-8].  The  electrochemical  impedance 
spectroscopy  (EIS)  method  has  been  used  to  characterize  the  bat¬ 
tery  impedance  behavior  over  a  wide  frequency  range  [5,9].  It  is 
also  recognized  to  have  the  ability  to  distinguish  different  physical 
processes  inside  batteries  and  determine  their  process  rates  by 
non-destructive  measurements  [10].  When  using  EIS,  a  properly 
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designed  battery  impedance  model  is  needed.  This  impedance 
model  is  expected  to  depict  battery  internal  dynamics,  including 
charge-transfer  reaction  on  electrode/electrolyte  interface,  lithium- 
ion  diffusion  in  electrodes,  double  layer  effects,  and  the  resistance/ 
capacitance  growth  of  the  anode  insulating  film,  etc.  [11]. 

Many  earlier  impedance  models  are  equivalent  circuit  typed 
because  of  their  simple  form  12].  The  most  classic  one  is  the 
Randles  model,  which  includes  a  Warburg  element,  a  resistor  and  a 
RC  circuit.  The  common  understanding  is  that  the  Warburg  element 
represents  diffusion  dynamics,  the  resistor  represents  pure  ohmic 
resistance,  while  the  RC  component  represents  the  double  layer 
effects.  The  Randles  model  and  its  variants  have  been  used  to  es¬ 
timate  state  of  health  (SOH),  e.g.  Waag  et  al.  [13],  Eddahech  et  al. 
[14]  and  Mukoyama  et  al.  7].  One  main  disadvantage  is  that  it  is  not 
easy  to  directly  relate  internal  physical  properties  of  batteries  to 
model  prediction.  This  disadvantage  has  limited  the  application  of 
this  type  of  models  on  designing  high-performance  battery 
management  systems  and  battery  designs. 
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Significant  work  has  been  published  on  electrochemistry-based 
impedance  expression.  Meyers  et  al.  first  proposed  an  impedance 
model  for  a  single  particle  and  then  extended  it  for  electrodes  with 
porous  structure  [15].  The  modeled  dynamics  include  the  charge- 
transfer  reaction,  double  layer  effect  and  solid-phase  diffusion. 
The  particle-size  distribution  is  considered  to  accurately  describe 
the  solid-phase  diffusion.  This  model  ignored  the  electrolyte 
diffusion  limitations  for  simplicity,  which  might  limit  its  applica¬ 
tion  to  high  current  charging/discharging  conditions.  In  addition,  it 
is  assumed  in  Meyers’s  study  that  there  is  only  one  phase  in  elec¬ 
trode,  however,  some  researchers  argued  that  certain  types  of 
lithium-ion  batteries  may  have  multiple  phases,  e.g.  Andersson 
et  al.  [16,17].  For  two-phase  particles,  Huang  et  al.  developed  a 
semi-mathematical  impedance  model  by  combining  electro¬ 
chemical  components  and  equivalent  circuit  components  18]. 
Huang’s  model  considers  two  homogenous  phases,  which  is  better 
suited  for  lithium  iron  phosphate  (LFP)  batteries.  Sikha  et  al.  further 
considered  the  diffusion  process  in  electrolyte  when  modeling  an 
insertion  cathode/separator/foil  anode  cell  [19].  This  work  is  also 
extended  to  a  dual  insertion  electrode  battery,  in  which  the  inser¬ 
tion  electrodes  are  separated  by  an  ionic  conductive  membrane 
[20].  Abraham  et  al.  developed  an  electrochemical  model  to  fit  the 
impedance  data  of  LiNio.sCoo.isAlo.osC^-based  positive  electrodes, 
in  which  the  diffusion  dynamics  in  electrolyte,  solid— electrolyte 
interface  (SEI)  and  active-material  surface  layer  was  considered 
[21  ].  These  electrochemistry-based  impedance  models  can  describe 
battery  internal  dynamics  such  as  ion  diffusion  and  migration  in 
electrodes  and  electrolyte,  charge-transfer  reaction  on  the  elec¬ 
trode/electrolyte  interface,  etc.  However,  the  relatively  complex 
structure  and  large  number  of  parameters  are  major  challenge  in 
practical  battery  management  systems.  As  far  as  we  know,  only  a 
few  models  have  been  used  on  some  narrow  topics,  e.g.,  how  to 
study  and  predict  the  aging  effect  [5].  In  addition,  direct  and 
quantitative  relationship  between  electrochemistry  and  impedance 
models  is  not  well  established  because  some  important  electro¬ 
chemical  features  have  been  neglected,  including  electrode 
porosity  and  the  variation  of  ion  concentration.  Bridging  this  gap 
will  help  to  understand  battery  dynamics  and  design  a  more  effi¬ 
cient  and  effective  energy  storage  system. 

The  main  purpose  of  this  paper  is  to  derive  a  structurally  concise 
but  accurate  battery  impedance  model  that  directly  relates  elec¬ 
trochemical  and  physical  parameters  with  battery  internal 
dynamics.  This  paper  focuses  on  lithium-ion  batteries,  beginning 
with  the  well-known  single  particle  assumption  followed  by  en¬ 
hancements  for  better  tradeoff  between  accuracy  and  simplicity. 
The  remainder  of  this  paper  is  structured  as  follows:  Section  2 
discusses  the  structural  simplification  for  lithium-ion  batteries.  In 
Section  3,  an  electrochemistry-based  impedance  model  is  derived. 
Section  4  describes  the  EIS  experiments  used  for  model  validation, 
followed  by  to  the  methodology  used  for  model  parameter  iden¬ 
tification.  In  Section  5,  the  battery  aging  effect  is  discussed  using 
identified  parameters  at  different  stage  of  aging.  Section  6  con¬ 
cludes  this  paper. 

2.  Simplified  battery  electrode  model 

This  paper  focuses  on  batteries  with  lithium  oxide  based  cath¬ 
ode  and  lithium-carbon  based  anode.  In  such  batteries,  the  charge- 
transfer  reaction  happens  on  the  electrode/electrolyte  interface. 
During  discharging,  the  lithium-ions  de-intercalate  from  anode, 
move  through  the  electrolyte,  and  intercalate  into  cathode.  When 
charging,  the  process  is  reversed.  To  tradeoff  between  accuracy  and 
conciseness,  it  is  assumed  that  (1)  electrodes  consist  of  multiple 
spherical  particles,  and  all  particles  have  identical  size  and  dy¬ 
namics,  (2)  the  material  inside  each  particle  is  spatially 


homogeneous,  (3)  electrodes  have  perfect  conductivity,  i.e.,  no 
voltage  drop  inside  a  particle  or  among  particles,  (4)  the  SEI  film 
growth  caused  by  aging  only  happens  on  the  anode  surface,  (5)  the 
electrolyte,  separator  and  current  collectors  are  treated  as  lumped 
pure  resistors.  With  these  assumptions  each  electrode  is  repre¬ 
sented  by  one  lumped  spherical  particle,  shown  in  Fig.  1. 

In  addition  to  diffusion  dynamics  inside  electrode,  charge- 
transfer  reaction  and  double  layer  effect,  which  are  frequently 
considered,  we  consider  the  following  behaviors  in  our  model. 

(1)  The  relationship  of  lithium-ion  concentration  and  the  open- 
circuit  potential:  the  lithium-ion  concentration  fluctuates 
when  charging  and  discharging,  which  leads  to  variations  in 
electrochemical  potential  and  overall  battery  open-circuit 
potential.  The  open-circuit  potential  is  often  empirically 
fitted  as  a  mathematical  function  of  SOC,  e.g.  linear  function 
[15],  higher-order  polynomials  [19,20  ,  etc.  These  functions 
do  not  relate  to  battery  dynamics  and  provide  no  insight  to 
electrochemical  dynamics.  We  choose  to  use  the  Nernst 
equation  directly  to  bridge  this  gap  from  the  viewpoint  of 
thermodynamics. 

(2)  The  capacitance  dispersion  caused  by  the  porosity  of  elec¬ 
trodes:  both  electrodes  at  present  are  porous  to  increase  the 
effective  surface  area.  This  porosity  leads  to  battery  imped¬ 
ance  much  more  complicated  than  any  commonly  used 
equivalent  circuit  models.  In  theory,  the  porous  electrodes 
can  be  approximated  by  multiple  spherical  particles  with 
predefined  size  distribution.  This  approximation,  however, 
often  causes  unacceptable  difficulty  in  model  structure  and 
parameter  identification.  Here,  a  constant  phase  element 
(CPE)  component  will  be  used  to  describe  the  capacitance 
dispersion  in  the  double  layer  effect  [9]. 

(3)  The  SEI  film  growth  at  anode  during  degradation:  the  SEI  film 
growth  is  a  major  mechanism  of  battery  degradation.  The 
film  growth  reduces  cyclable  lithium-ions  and  increases 
internal  resistance,  which  leads  to  both  capacity  loss  and 
power  fade  which  is  an  irreversible  side  reaction,  especially 
on  the  anode  side  [22,23].  This  paper  ignores  the  SEI  film 
growth  on  the  cathode  side  because  of  its  lower  significance 
[24]. 

Fig.  2  sketches  the  current  flow  during  discharging  for  both 
anode  and  cathode.  The  current  flow  at  cathode  is  similar  to  that  of 
anode,  just  without  the  components  related  to  the  SEI  film. 

In  Fig.  2,  Jfd.i  represents  the  Faradaic  current  density  due  to 
charge-transfer  reactions,  Zfd.i  is  the  Faradaic  impedance  of 
electrodes  caused  by  charge-transfer  reactions  (i  =  n  for  anode  and 
i  =  p  for  cathode),  determined  from  the  Faradaic  current  density 
and  the  potential  across  the  interface,  0 y,  02, i  represent  the  po¬ 
tential  on  the  two  sides  of  electrode/electrolyte  interface,  Jdi.i  reP_ 
resents  the  double  layer  current  density,  Cdi,/  represents  the  double 
layer  capacitance,  Cmm  represents  the  capacitance  due  to  the  SEI 


Anode 


Separator  Cathode 


Anode  Separator  Cathode 


Fig.  1.  Simplification  of  lithium-ion  battery  electrodes  (left  -  multiple  particle;  right  - 
single  particle). 
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Fig.  2.  Current  flows  at  the  anode  and  cathode  ((a)  anode;  and  (b)  cathode). 


film  insulation  on  anode  Jmm  is  the  current  density  through  the  SEI 
film,  and  Rmm  is  the  SEI  film  resistance.  In  general,  the  impedance 
spectrum  of  lithium-ion  batteries  includes  multiple  sections,  for 
instance,  a  straight  tail  at  low  frequency,  a  semi-circle  at  middle 
frequency  and  a  quarter-circle  at  high  frequency.  Detailed  discus¬ 
sion  will  be  in  Section  4. 

3.  Battery  impedance  model 

This  section  first  reviews  the  partial  differential  equations 
(PDEs)  used  to  depict  battery  electrochemical  dynamics,  including 
charge-transfer  dynamics,  open-circuit  potential  and  ion  diffusion 
dynamics  in  the  electrodes.  Then,  we  derive  the  Faradaic  imped¬ 
ance  from  the  well-known  Butler-Volmer  equation,  Nernst  equa¬ 
tion  and  Fields  second  law,  followed  by  impedance  models  for 
anode  and  cathode.  The  overall  battery  impedance  is  obtained  by 
adding  other  components,  including  the  double  layer  capacitance, 
the  SEI  film  impedance  and  the  pure  inductance. 


3.1.  PDEs  for  the  single  particle  model 


3.1.1.  Electrochemical  dynamics  on  electrode  surface 

The  Faradaic  current  density  is  governed  by  the  well-known 
Butler-Volmer  equation,  which  describes  the  relationship  be¬ 
tween  the  overpotential  and  Faradaic  current  density,  expressed  as 


Jfd,i 


hi 


-  exp 


(1  -  «)F 

RT 


i  =  n,p, 


(1) 


Vi  =  hi  -<t>2,i  -Ei,  i  =  n,p,  (2) 

where  £z  is  the  open-circuit  potential  of  anode  and  cathode.  Here, 
we  use  an  expression  derived  from  the  Nernst  equation  instead  of 
linear  functions  or  higher-order  polynomials,  which  is  more  suit¬ 
able  to  describe  battery  inherent  mechanism  while  maintaining  a 
concise  structure  [26] 

L  =  Eo,i~Y ln(1  y  i  =  n,p,  (3) 

where  Fo,i  is  the  reference  open-circuit  potential,  Xi  is  the  mole 
fraction  of  the  lithium-intercalated  sites  in  the  electrodes.  The  mole 
fraction  is  defined  as 


Xi  =  ViC^i/cii,  i  =  n,p,  (4) 

where  C\  z  is  the  average  lithium-ion  concentration  of  electrodes,  Vj 
is  the  electrode  volume,  and  az  is  the  maximum  mole  number  of 
lithium-ions  in  the  electrodes.  Vz  is  often  assumed  to  be  constant 
when  a  relatively  small  alternative  perturbation  is  applied  in  EIS 
measurement  [27].  Let 

Iq  =  Vj/cii,  i  =  n,p,  (5) 


and  substitute  Eqs.  (4)  and  (5)  into  Eq.  (3),  we  have 


kiC 


r  c  RT, 

Ei  =  Eo,i--p-ln-  ,  _  , 

F  1  -  k,Cu 


i  =  n,p. 


(6) 


In  Eq.  (6),  further  information  is  required  to  determine  the 
average  lithium-ion  concentration  of  electrodes.  The  average 
lithium-ion  concentration  inside  a  particle  is  almost  equal  to  that 
on  the  surface  when  battery  current  is  close  to  zero  [28  .  Therefore, 
a  relationship  between  electrode  overpotential  and  bulk  ion  con¬ 
centration  can  be  built.  In  EIS  tests,  the  current  applied  is  relatively 
small  and  a  similar  assumption  can  be  used  in  the  following  model 
derivation. 


3.1.2.  Diffusion  dynamics  at  anode  and  cathode 

The  lithium-ion  concentration  in  both  electrode  and  electrolyte 
is  affected  by  the  ion  diffusion  dynamics.  When  modeling,  the 
diffusion  in  electrolyte  is  much  faster  than  that  in  electrode  and 
thus  is  neglected.  The  diffusion  of  lithium-ions  in  electrodes  is 
governed  by  the  Fields  second  law  [29] 


dCi ,j  _  9  /  2^1, A 

at  -  rf  9rz  \  1  9rz  J  ’ 


i  =  n,p, 


(7) 


where  rz  is  the  radial  position  in  the  spherical  coordinate,  Cy  de¬ 
notes  lithium-ion  concentration  in  the  particles  at  radius  rz  and 
time  t,  Di)Z  is  the  diffusion  coefficient.  The  boundary  conditions  are 

"h  =0,  i  =  n,p ,  (8) 


where  p z  is  the  overpotential,  a  is  the  transfer  coefficient  (a  =  0.5  as 
given  in  Ref.  25]),  Jo.z  is  the  exchange  current  density  in  the  nom¬ 
inal  condition,  F  is  the  Faraday’s  constant,  R  is  the  universal  gas 
constant  and  T  is  the  absolute  temperature.  In  normal  charging  and 
discharging  conditions,  Jo,i  is  not  a  constant,  but  depends  on  the 
lithium-ion  concentration  in  the  electrodes  and  electrolyte.  In  the 
EIS  analysis,  the  battery  is  excited  with  a  low  amplitude  current  and 
therefore  J0,i  is  approximately  constant  [15].  The  overpotential  rji  is 
defined  as 


dCu 

9rz 


Ti=Ri 


J  FD,z 
DUF’ 


i  =  n,p. 


(9) 


The  boundary  conditions  of  Eq.  (8)  are  because  spherical  parti¬ 
cles  are  radially  symmetric  and  Eq.  (9)  is  from  the  mass  and  charge 
conversion  at  the  solid-electrolyte  interface.  It  is  reasonable  to 
ignore  the  electrolyte  diffusion  dynamics  when  the  charge  and 
discharge  rate  is  low  because  of  the  higher  diffusivity  in  the  elec¬ 
trolyte  [25,30].  This  condition  holds  true  for  our  EIS  tests. 


12 


S.E.  Li  et  al.  /  Journal  of  Power  Sources  258  (2014)  9-18 


3.2.  Derivation  of  Faradaic  impedance 


The  diffusion  dynamics  in  the  electrodes  are  described  by  partial 
differential  equations  (PDEs).  The  PDEs  are  rewritten  into  ordinary 
differential  equations  (ODEs)  using  the  Laplace  transformation 
method.  In  the  EIS  tests,  the  battery  is  excited  by  a  small  sinusoidal 
signal  and  all  battery  dynamics  works  in  linear  regions  around 
equilibrium  points.  Hence,  it  is  assumed  that  both  the  changes  of 
lithium-ion  concentration  and  Faradaic  current  density  fluctuates 
periodically  at  the  same  frequency 


JfD,i  -  JfD,z  +  JfD,!6^  I  =  v  n 

Cu(r,t)  *£U  +  Cu(r)ei'*' 


(10) 


where  variables  with  a  “bar”  represents  the  nominal  values  and 
variables  with  a  “tilde”  represent  the  amplitude  of  perturbation,  w 
is  the  excitation  frequency.  Plugging  Eq.  (10)  to  Eq.  (7),  we  have 


juCu 


Phi  Q  />Ci,A 
rf  9 n  (  '  ar,  )' 


i  =  p,n. 


(11) 


Several  parameters  in  Eq.  (17),  namely,  CitStj,  Ru  D y  and  /q  are 
coupled  and  cannot  be  identified  separately.  Hence,  we  define  two 
lumped  parameters 

^l.i  —  Ci)S,A/(l  —  ^iQ,s,i)/^L  *  = 

K2,i  =  Ri/(VDv)  1  J 

Then  the  Faradaic  impedance  of  unit  surface  area  is  then 
expressed  as 

z  bfay  Kt(  1  _ 1 _ | 

AjFDi  F  V01 

i  =  n,p. 

(19) 

The  Faradaic  impedance  has  two  items,  charge-transfer  resis¬ 
tance  and  the  impedance  caused  by  solid-phase  diffusion.  The 
charge-transfer  resistance  is  a  critical  parameter  in  the  battery 
behavioral  analysis.  Let  Rct,i  denote  the  charge-transfer  resistance  of 
unit  surface  area,  expressed  as 


Now  Eq.  ( 11 )  becomes  a  nonlinear  second  order  ODE.  Solving  Eq. 
(11)  together  with  the  associated  boundary  conditions  Eqs.  (8)  and 
(9),  we  obtain  the  lithium-ion  concentration  on  the  surface  of 
particles  (r;  =  R{) 


fl  ,s,i  —  Ct,i(Ri)  —  01,JfD,i 


h.i  = 


i  =  n,p. 


(12) 


In  the  following,  we  use  Ci>s,i  to  denote  the  lithium-ion  con¬ 
centration  on  the  surface  of  particles.  Substituting  Eq.  (6)  into  the 
linearized  Butler— Volmer  equation,  we  have 


Jfdj 


1  -  frfCljSjj )  ’ 


i  =  n,p, 


(13) 


where  (pu,i  =  0i,z  -  02, u  Differentiating  Eq.  (13), 


4/fd,z  =  hixfWnj  3- Joy 


’'Ci,s,i(l  -  fc*clA,) 


ACljSji,  i  =  n,p. 


(14) 


When  excited  by  a  small  sinusoidal  signal,  the  variations  of  all 
variables  in  Eq.  (14)  are  also  sinusoidal.  Hence,  we  have 


Rcti  =  RT/(joiF),  i  =  p,n.  (20) 

The  overall  Faradaic  impedance  is  related  to  the  surface  area. 
The  following  model  identification  process  will  identify  the  whole 
Faradaic  impedance,  rather  than  that  per  unit  surface  area. 

3.3.  Derivation  of  electrode  impedance 

Besides  the  Faradaic  impedance,  the  electrode  impedance  in¬ 
cludes  the  effects  of  double  layer  and  SEI  film  (the  cathode  SEI  film 
impedance  is  assumed  to  be  zero  in  this  paper).  We  will  lump  them 
together  and  derive  the  impedance  model  for  cathode  and  anode. 

3.3 A.  Cathode  impedance  per  unit  area 

Fig.  2(b)  shows  the  current  flows  in  the  cathode.  The  schematic 
diagram  of  its  equivalent  circuit  structure  is  shown  in  Fig.  3. 

The  porous  electrodes  have  rough  electrode  surface  which  lead 
to  double  layer  effect  with  capacitance  dispersion,  rather  than 
behaving  like  an  ideal  capacitor  [31].  A  constant  phase  element 
(CPE)  is  used  to  describe  this  dispersed  capacitance  [9],  expressed  as 

zdu  *  l/(cOua<'),  i  =  P, n,  (21) 

where  Zdi,/  is  the  impedance  of  the  CPE,  u*  is  the  exponent  to 
describe  the  dispersion  effect  (0  <  v  *<  1 ).  When  u*  =  1,  CPE  becomes 
an  ideal  capacitor.  Combined  with  Eq.  (21),  the  cathode  impedance 
of  unit  area  is 


4/fD,z  —  JFD,i  -JfDJ  —  JFD,ieJ0)t 

ACi,Sii  =  =  ClA,efr*  i  =  n,p.  (15) 

A(/>i2  =  <t>  12  —  ^12  =  l^12e,Wt 

Substitute  Eqs.  (12)  and  (15)  into  Eq.  (14),  we  have 


J  FD,z 


Jo.iF  ^  ,  Jo.i 

W  012  C1)S)l-(l-lqCi)S)l-) 


0i,Jfd,i,  i  =  n,p. 


(16) 


The  Faradaic  impedance  of  unit  surface  area  is  then  expressed  as 


7  _  Nhij  _  RT 

ZFD  ,i  -  aj^7  _  ~r 


Ci,S)i(l 


FD, 


i  =  p,  n 


(17) 


2p  =  Zfd,p/(i  +  ^FD,pCd|  p(Jw)Up)  •  (22) 

3.3.2.  Anode  impedance  of  unit  area 

The  anode  particle  has  more  complicated  equivalent  structure 
as  shown  in  Fig.  4.  The  CPE  is  used  for  both  double  layer  capaci¬ 
tance  and  SEI  film.  Here  we  assume  that  the  SEI  film  grows 
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uniformly  over  the  porous  electrode,  as  a  result  of  which,  the  CPE 
representing  SEI  film  impedance  has  the  same  exponent  as  anode 
double  layer  capacitance.  The  anode  impedance  of  unit  area,  Zn  sEi 
is  expressed  as 


Zn,SEI 


_ Zn  +  flfilm _ 

1  +  0'<nCglm(Zn  +  Rfilm)’ 


(23) 


where  Rf\\m  is  the  resistance  of  SEI  film,  and  Cj?lm  is  a  capacitance¬ 
like  parameter. 


3.4.  Full  battery  impedance  model 

Fig.  5  shows  the  model  structure  of  the  whole  battery.  The 
model  structure  has  an  anode  impedance  and  a  cathode  impedance 
connected  in  series.  The  impedance  of  electrolyte,  separator,  and 
current  collector  are  represented  by  a  pure  ohmic  resistor  %  placed 
between  anode  and  cathode.  At  high  frequency,  battery  current 
collectors  and  cables  introduce  inductance,  which  appears  as  a 
positive  imaginary  part  in  the  EIS  diagram.  This  inductive  effect  is 
represented  by  a  pure  inductor  L  in  Fig.  5. 

The  overall  impedance  of  a  battery  Zbattery  then  has  the  following 
form: 

^battery  —  jwL  +  ZnSEl/Sn  +  R-o  +  Zp/Sp,  (24) 

where  Sp  and  Sn  are  the  cathode  and  anode  effective  surface  area, 
respectively. 


the  impedance  curve  gradually  shifts  rightward  and  the  semi¬ 
circle  at  middle  frequency  enlarges  with  the  increasing  cycle 
number,  and  the  enlarging  semi-circle  swallows  the  straight  tail  at 
low  frequency.  This  is  an  obvious  index  to  show  the  increase  of 
Faradaic  impedance  and  the  slowdown  of  interfacial  charge- 
transfer  reaction.  The  state  of  charge  also  affects  the  shape  and 
amplitude  of  impedance.  It  is  observed  that  the  impedance  of  half 
charged  batteries  is  smaller  than  that  of  the  fully  charged 
batteries. 

4.2.  Model  parameter  identification 

The  above-derived  impedance  model  has  16  parameters,  listed 
in  Tables  2  and  3,  respectively.  These  parameters  are  closely  related 
to  battery  physical  and  electrochemical  characteristics.  An  identi¬ 
fication  method  is  used  to  minimize  the  error  between  impedance 
measurement  and  model  prediction.  The  cost  function  to  be  opti¬ 
mized  is 


4.  Impedance  measurement  and  parameter  identification 

4.1.  Experiment  setup  and  EIS  tests 

The  EIS  measurement  was  conducted  in  a  battery  impedance 
tester,  which  has  two  channels,  capable  of  measuring  two  cells 
simultaneously.  The  two  battery  cells  are  lithium-ion  typed.  Its 
specifications  are  listed  in  Table  1. 

The  impedance  tests  are  carried  out  for  two  cells  (noted  as  A 
and  B)  at  ambient  temperature  (25  °C)  and  two  SOC  levels  (half 
charged  SOC  ~50%  and  fully  charged  SOC  ~100%).  Fig.  6  shows 
the  impedance  measurement  of  both  Cell  A  and  Cell  B.  The  solid 
lines  represent  different  aging  cycles,  200,  400,  600,  800  and 
900. 

The  measured  impedance  contains  three  parts:  a  straight  tail  at 
low  frequency,  a  large  semi-circle  at  middle  frequency,  and  a  small 
quarter-circle  at  high  frequency.  The  general  understanding  is  that 
the  straight  tail  relates  to  the  diffusion  dynamics  of  lithium-ions  in 
electrodes,  the  larger  semi-circle  at  middle  frequency  is  associated 
with  the  interfacial  charge-transfer  reaction  combined  with  double 
layer  capacitance,  and  the  quarter-circle  at  high  frequency  is  due  to 
the  presence  of  SEI  film  and  the  particle-to-particle  contact  of 
active  materials  [13]. 

The  aging  cycles  have  a  direct  influence  on  the  shape  and 
amplitude  of  impedance.  It  is  observed  that  as  battery  degrades, 


/max  .  /  ^\\  2 

min  e=  £  (zexp(f)  -ZhMeJf,d))  ,  (25) 

0  /=/min 

where  6  is  the  estimated  parameter  set,  e  is  the  model  error,  Zex p  is 
the  measured  data  and/is  the  frequency  point.  The  Hybrid  Multiple 
Particle  Swarm  Optimization  method  (HMPSO)  [32]  is  used  to 
obtain  the  optimal  parameters.  The  HMPSO  is  a  global  optimization 
method  with  effective  memory  utilization,  more  robust  for  noisy 
data  and  more  efficient  in  maintaining  diverse  set  of  search  can¬ 
didates.  The  effectiveness  has  been  demonstrated  in  battery 
equivalent  circuit  model  studies  [33]. 

Before  running  the  optimal  parameter  identification,  we 
reduced  the  over-fitting  problem  by  fixing  some  parameters  which 
are  less  sensitive  to  battery  SOC  and  aging.  The  fixed  parameters  are 
listed  in  Table  2.  The  reasons  to  fix  these  6  parameters  are  explained 
below.  According  to  [5],  Cj  ,  Cjn  and  Sp,  Sn  are  not  sensitive  to 
model  error  and  are  regarded  as  constant.  We  also  decide  to  fix  I<^n 
and  L  based  on  the  observation  of  the  identification  results.  The 
sensitivity  analysis  by  varying  these  two  parameters  showed  less 
effect  on  model  fitness  than  other  parameters,  moreover,  fixing 
these  two  parameters  helps  to  avoid  over-fitting.  Such  difficult 
parameters  have  been  explained  by  Troltzsch  [5].  The  values  of  Sp 
and  Sn  are  selected  based  on  known  average  values  of  electrode 
surface  areas.  The  values  of  C°lp,  Cjjln  /< \n  and  L  are  determined 
from  results  of  all  aging  cycles. 


Table  1 

Specifications  of  lithium-ion  battery  cells. 


Variable 

Value 

Capacity 

940  mAh 

Nominal  voltage 

3.7  V 

Max.  voltage 

4.2  V 

Min.  voltage 

2.5  V 

Energy  density 

133  Whkg"1 

Fig.  4.  Anode  equivalent  circuit  structure. 
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Fig.  6.  Measured  impedance  under  different  SOC  and  aging  cycles  ((a)  Cell  A  -  fully  charged;  (b)  Cell  A  -  half  charged;  (c)  Cell  B  -  fully  charged;  and  (d)  Cell  B  -  half  charged). 


Table  2 

Fixed  parameters. 


No. 

Parameters 

Description 

Value 

Unit 

1 

2 

C° 

> 

'“dl,n 

Cathode  double  layer  capacitance 
Anode  double  layer  capacitance 

7.88 

1.98 

F  sup-1  m-2 

F  sUn_  1  m-2 

3 

sp 

Cathode  effective  surface  area 

3.52 

m2 

4 

Sn 

Anode  effective  surface  area 

2.16 

m2 

5 

K l,n 

Lumped  parameter  of  anode 

1.63 

Am-2 

6 

I 

Pure  inductor 

1.7e-07 

H 

4.3. 

Model  accuracy  analysis 

The  Randles  model  is  selected  to  be  the  baseline  for  comparison. 
Figs.  7  and  8  illustrate  the  fitting  results  of  two  models  (Randles 
model  and  our  proposed  electrochemistry-based  impedance 
model).  Figs.  9  and  10  compare  their  accuracy  based  on  root  mean 
square  error  (RMSE)  using  both  Cell  A  and  Cell  B  data. 

Here,  we  will  discuss  where  the  Randles  model  loses  accuracy 
and  how  an  electrochemistry-based  model  makes  improvement.  It 
is  observed  from  the  identification  results  that  the  Randles  model 


Table  3 

Free  parameters. 


No. 

Parameters 

Description 

Unit 

7,8 

Jo,u  i  =  P,n, 

Exchange  current  density 

A  m2 

9 

Ki  ,P 

Lumped  parameter  of  cathode 

A  m2 

10,11 

K 2,i,  i  =  p,n, 

Lumped  parameter  of  electrodes 

s0-5 

12,13 

vu  i  =  p,n, 

Exponent  of  CPE  element 

- 

14 

Kfilm 

SEI  film  resistance 

O  m2 

15 

C° 

'“film 

SEI  film  capacitance 

F  sUn_1  m-2 

16 

Ro 

Pure  resistance 

Q 

loses  accuracy  in  all  three  frequency  range.  In  the  low  frequency, 
the  Randles  model  is  bound  to  have  a  45  degree  tail,  but  actual  EIS 
data  have  a  slope  higher  than  45  degree.  We  think  this  is  because 
ion  diffusion  in  the  electrodes  is  in  spherical  structure  and  is  not 
along  a  single  dimension.  The  latter  is  the  assumption  of  Warburg 
element  used  in  the  Randles  model.  In  the  middle  frequency,  the 
Randles  model  predicts  the  impedance  to  be  a  semi-circle,  which  is 
different  from  what  actual  EIS  data  look  like.  The  main  reason  is  the 
capacitance  dispersion  due  to  the  porosity  of  electrodes.  In  the 
electrochemistry-based  model,  the  CPE  component  can  describe 
this  frequency-dependent  characteristic,  which  demonstrates  the 
importance  of  considering  electrode  porosity.  The  Randles  model 
loses  the  ability  to  depict  high  frequency  impedance  because  it 
omits  some  impedance  components.  The  model  accuracy  on 
average  is  0.65  mQ  (RMSE,  root  mean  square  error)  for  the  pro¬ 
posed  model,  about  one-eighth  that  of  the  Randles  model. 

5.  Parametric  variation  due  to  aging 

It  is  apparent  that  battery  impedance  changes  significantly  with 
aging.  The  impedance  curve  gradually  shifts  rightward  and  the 
semi-circle  at  the  middle  frequency  region  enlarges  with  the  bat¬ 
tery  age.  We  identified  the  battery  parameters  of  the  proposed 
electrochemistry-based  model  with  the  abovementioned  EIS  data. 
The  identified  parameters  of  two  cells  are  averaged  at  the  same  SOC 
and  aging  cycle  in  order  to  reduce  the  uncertainties  of  cell  incon¬ 
sistency.  Table  4  lists  the  average  critical  parameters  of  fully 
charged  state  (at  different  aging  cycles:  200,  400,  600,  800  and 
900),  while  Table  5  lists  those  of  half  charged  state  (at  different 
aging  cycles:  200,  400,  600,  800  and  900). 

Fig.  11  shows  the  relationship  between  pure  ohmic  resistance 
and  aging  cycles.  Some  researchers  believed  that  the  pure  ohmic 
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Re(Z)(Q) 


Re(Z)(Q) 


Re(Z)(Q)  Re(Z)(Q) 


Fig.  7.  Fitting  results  of  battery  impedance  (Cell  A)  ((a)  Cycle  200,  fully  charged;  (b)  Cycle  200,  half  charged;  (c)  Cycle  800,  fully  charged;  and(d)  Cycle  800,  half  charged). 


Re(Z)(Q) 


Re(Z)(Q) 


Re(Z)(Q) 


Re(Z)(Q) 


Fig.  8.  Fitting  results  of  battery  impedance  (Cell  B)  ((a)  Cycle  200,  fully  charged;  (b)  Cycle  200,  half  charged;  (c)  Cycle  800,  fully  charged;  and  (d)  Cycle  800,  half  charged). 
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Fig.  9.  Comparison  of  model  accuracy  based  on  RMSE  (Cell  A)  ((a)  proposed  model;  and  (b)  Randles  model). 


Cycle  number 


Fig.  10.  Comparison  of  model  accuracy  based  on  RMSE 

resistance  Ro  is  equal  to  the  value  at  the  intersection  of  the  EIS 
measurement  with  the  real  axis  [6,22].  Waag  et  al.  argued  that  the 
intersection  with  the  real  axis  contains  pure  resistors  other  than  Ro 
[13].  This  phenomenon  is  supported  by  our  identification  results, 
i.e.  the  identified  Ro  is  smaller  than  the  corresponding  intersecting 
value,  by  comparing  Table  4  and  Fig.  6(c).  With  aging  from  cycle  200 
to  900,  R0  of  fully  charged  batteries  increases  from  62.7  mQ  to 
71.8  mQ,  while  that  of  half  charged  batteries  increases  from 
61.1  mQ  to  73.9  mQ.  This  is  mainly  caused  by  the  decomposition  of 
electrolyte,  as  described  in  Ref.  [22].  In  addition,  the  pure  ohmic 
resistance  of  fully  charged  and  half  charged  batteries  is  almost 
identical  at  the  same  aging  level,  implying  that  the  pure  ohmic 
resistance  does  not  vary  with  SOC,  which  is  also  concluded  in  other 
publications  [13]. 

As  seen  in  Fig.  12,  the  SEI  film  resistance  Rmm/Sn  increases  with 
the  cycle  number.  For  fully  charged  state,  Rmm/Sn  increases  from 
2.3  mQ  to  6  mCl  for  aging  cycle  200-900,  almost  tripled  in  ampli¬ 
tude.  Meanwhile,  Rmm/Sn  increases  from  2.8  mQ  to  5.2  mQ  for  half 


Cycle  number 


(Cell  B)  ((a)  proposed  model;  and  (b)  Randles  model). 

charged  cells  for  aging  cycle  200-900,  almost  doubled  in  ampli¬ 
tude.  It  is  believed  that  the  growth  of  anode  SEI  is  a  key  aging 
mechanism,  leading  to  battery  power  and  capacity  degradation 
[6,22,34-36].  The  film  resistance  growth  can  be  a  measure  for 
predicting  capacity  and  power  fade  of  lithium-ion  batteries.  The 
increasing  trend  of  Rmm/Sn  is  not  as  consistent  as  that  of  R0. 
Moreover,  its  value  is  affected  by  battery  SOC  levels.  The  time 
constant  imm  is  the  product  of  Rmm  and  C{~lm.  As  shown  in  Fig.  13, 
the  parameter  imm  also  increases  with  battery  age.  It  is  observed 
that  Tfiim  at  half  charged  state  is  larger  than  that  at  fully  charged 
state. 

In  a  lithium-ion  battery,  the  charge-transfer  resistance  in  par¬ 
allel  with  the  double  layer  capacitance  dominates  the  depressed 
semi-ellipse  at  the  mid  frequency.  In  a  two-electrode  system,  it  is 
hard  to  distinguish  the  charge-transfer  resistances  of  anode  and 
cathode.  Here,  the  charge-transfer  resistance  on  the  cathode  side  is 
fixed  in  order  to  better  demonstrate  the  aging  effects  at  the  anode 
side.  When  aging  from  cycle  200-900,  the  value  of  Rct,n/5n( fully 


Table  4 

Identified  parameters  vs.  cycle  number  (average  of  fully  charged  state). 


Cycle 

Number 

Ro  (O) 

Rct,n/S„(Q) 

WSn(O) 

'Oilm  (s) 

K2,n  (s0-5) 

200 

6.27e-02 

1.21e-02 

2.34e-03 

5.38e-04 

3.19 

400 

6.59e-02 

3.04e-02 

3.84e-03 

9.23e-04 

3.13 

600 

6.65e-02 

6.46e-02 

6.12e-03 

1.12e-03 

3.01 

800 

7.18e-02 

6.09e-02 

6.33e-03 

2.74e-03 

3.64 

900 

7.18e-02 

7.54e-02 

5.99e-03 

2.50e-03 

3.26 

Table  5 

Identified  parameters  vs.  cycle  number  (average  of  half  charged  state). 


Cycle 

Number 

Ro  (O) 

Rct,n/S„(Q) 

WS„(Q) 

^film  (s) 

I<2 ,n  (S°-5) 

200 

6.11e-02 

3.20e-03 

2.80e-03 

2.25e-03 

2.45 

400 

6.52e-02 

1.14e-02 

2.94e-03 

1.92e-03 

3.00 

600 

6.77e-02 

3.33e-02 

5.00e-03 

3.51e-03 

2.74 

800 

6.89e-02 

3.48e-02 

4.73e-03 

4.97e-03 

2.86 

900 

7.39e-02 

4.34e-02 

5.23e-03 

4.94e-03 

2.94 
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Fig.  12.  SEI  resistance  Kfiim/Sn  vs.  cycle  number. 


The  parameter  Tw,n,  is  defined  to  be  equal  to  the  square  of 
lumped  parameter  I<2, n  (where  n  denotes  anode): 

Tw,n  =  (I<2,  n)2  =  fi/Dj, (26) 

From  the  viewpoint  of  thermodynamics,  Tw>n  represents  the 
time  constant  of  the  diffusion  dynamics  on  the  anode  side.  Fig.  15 
depicts  the  relationship  between  Tw,n  and  the  aging  cycles.  It  is 
found  that  TW)n  for  both  half  and  fully  charged  batteries  show  no 
obvious  increasing  or  decreasing  tendency  with  cycle  number. 
Moreover,  Tw,n  of  half  charged  batteries  is  smaller  than  that  of  fully 
charged  battery.  A  common  understanding  is  that  the  particle 
radius  keeps  constant  at  different  aging  and  SOC  levels.  Any  change 
in  iWin  can  be  related  to  that  of  diffusion  coefficient  Di  n.  Therefore, 
it  can  be  inferred  that  the  diffusion  coefficient  of  half  charged 
batteries  is  higher  than  that  of  fully  charged  batteries,  and  battery 
degradation  does  not  slow  down  the  diffusion  dynamics  of  lithium- 
ions.  Similar  summary  can  be  found  from  Ref.  [37]  for  cathode  and 
Ref.  [38]  for  anode. 


charged)  increases  from  12.1  mQ  to  75.4  mQ,  almost  by  6  times.  For 
half  charged  state,  the  value  of  Rct,n/Sn  increases  from  3.2  mQ  to 
43.4  mQ.  One  explanation  for  this  obvious  increase  is  that  the  SEI 
film  growth  reduces  charge-transfer  reaction  rate  at  the  electrode/ 
electrolyte  interface.  In  addition,  it  is  found  from  Fig.  14  that  the 
charge-transfer  resistance  at  half  charged  state  is  smaller  than  that 
of  fully  charged  state.  This  phenomenon  can  be  explained  by  the 
change  of  exchange  current  density  caused  by  the  solid-phase 
lithium-ion  concentrations  [29]. 


6.  Conclusion 

An  electrochemistry-based  impedance  model  is  developed  to 
capture  the  impedance  behavior  of  lithium-ion  cells.  The  modeled 
dynamics  include  the  charge-transfer  reaction  at  the  electrode 
surface,  the  diffusion  dynamics  in  the  electrodes,  the  double  layer 
effects,  and  the  resistance/capacitance  growth  caused  by  anode 
insulating  film  in  degradation.  In  addition,  the  Nernst  equation  is 
directly  employed  to  describe  the  lithium-ion  concentration’s  effect 
on  open-circuit  voltage  and  a  constant  phase  element  is  used  to 
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Fig.  13.  Time  constant  Tfiim  vs.  cycle  number. 
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Fig.  15.  Time  constant  of  anode  diffusion  dynamics  vs.  cycle  number. 
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better  describe  the  capacitance  dispersion  caused  by  electrode 
porosity.  The  impedance  test  results  from  two  lithium-ion  battery 
cells  are  used  to  validate  the  proposed  model.  The  fitting  results 
show  that  this  model  works  well  throughout  the  aging  cycles.  Some 
identified  parameters  show  clear  trend  as  the  battery  ages,  which 
shows  good  potential  as  an  indicator  for  capacity  degradation.  The 
charge-transfer  resistance  and  the  SEI  film  resistance  demonstrate 
more  significant  trend  with  aging  cycles  and  the  results  are 
consistent  with  experiment-based  observations  from  the  literature, 
which  seems  to  be  indicating  the  potential  of  the  proposed  model 
for  battery  age  estimation. 
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